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O ■ I- INTRODUCTION AND MOTIVATION 

(N 

• In thepast year there have been many new interesting experimental discoveries in meson spectroscopy with heavy 
1^ ' quarks y . One of the goals of the experimental program at CLEO-c is to refine our experimental knowledge of 
l~H . the heavy hadron spectrum. This new data helps to validate methods, such as lattice gauge theory, that solve QCD 
^— I [ non-perturbatively. 

On the theoretical side, it has been claimed that there has been much progress in unquenched lattice QCD cal- 
culations [3- These new lattice QCD calculations that use improved staggered quarks have passed some important 
experimental consistency checks 3]. However, the one place where the agreement between lattice QCD and and 
' experiment is still poor is the mass splitting between the J ftp and rjc Q . The masses of these two mesons can usu- 
I ally be computed with the smallest statistical errors. Also, as these masses are independent of light valence quarks, 
■ this splitting does not depend on a large extrapolation in the valence quark mass. It does of course depend on an 
CnJ I extrapolation in the sea quark mass. 

. The experimental value for the mass splitting between the J/tp and ijc is 116 MeV. Some of the older lattice 
' calculations 1, 0, S H S ITol Ull that computed this splitting have been reviewed recently 531 • 

When using a clover improved fermion formalism on the lattice, the hyperfine splitting is sensitive to the coefficient 
csw (of ^ term in the fermion operator that helps reduce lattice spacing dependence on physical quantities) at nonzero 
lattice spacing, but the hyperfine splitting should be independent of the csw as the continuum limit is taken, because 
the clover term is an irrelevant operator. Recently the QCD-TARO collaboration have studied the charmonium 
' spectrum in quenched QCD using the clover action at a smaller lattice spacing {a~^ ^ 5 GeV) than previously used. 
They obtained a hyperfine splitting of 77(2) (6) MeV in the continuum limit. 
^ There has not been much work on the charmonium spectrum from unquenched lattice QCD calculations. There are 

• i-H I arguments based onpotential models, that suggest the hyperfine splitting in charmonium is sensitive to the presence 
' of light sea quarks [a, 01- El-Khadra et al. 17] did look at the charmonium spectrum on (unimproved) staggered 
^ ', gauge configurations {rriTr/mp was 0.6 and the lattice spacing was ^ 0.99(4) GeV.) from the MILC collaboration. 
No significant increase in the hyperfine splitting was reported. Stewart and Koniuk |18| studied the charmonium 
spectrum using NRQCD on unquenched (unimproved) staggered gauge configurations (m,r/wp ~ 0.45 and a ~ 0.16 
fm). Any signal for the effect of unquenching was hidden beneath the systematic uncertainties in using the NRQCD 
formalism for charmonium. 

The work by Davies et al. 0,^3 found that the correct ratio was produced for the (P-S)/(2S-1S) mass splittings for 
bb using unquenched calculations with improved staggered fermions. However, the hyperfine splitting in charmonium 
is still incorrect 0, l20l l2l| at 97 ± 2 MeV. The authors claim that the discrepancy may be caused by the clover 
coefficient only being used to tree level in tadpole imp roved perturbation theory. These calculations do, however, get 



the hyperfine splitting, the so-called J parameter |23|, in the light quark sector correct |19|. 

There is another possible reason that the hyperfine mass splitting between the J/ 4> and ri r is smaller than exper- 
iment in current lattice evaluations. This has been also discussed recently by QCD-TARO l2^ |2^. All lattice 
calculations have computed the non-singiet correlator (see figure [S]). However, charmonium interpolating operators 
are actually singlet (cFc), so the Wick contractions contain bubble diagrams (see figure (SJ. The bubble diagrams are 
OZI suppressed so should be small. However, this argument will fail if there is additional non-perturbative physics. 
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Group 


Method 


Md* - Mr, MeV 


Boyle f2_6] 
Boyle [22 
Hein et al. [28j 
UKQCD [29] 
UKQCD [30] 


clover 
(3=6.0 tadpole clover 
NRQCD P = 5.7 
NP clover /3 = 6.2 
NP clover /3 = 6.2 


124(8)(15) 
106(8) 
110lgtg^(3)(6)(5) 

127(14)(1)(3) 


PDG [25j 


Experiment 


142.12(7) 



TABLE I: Collection of hyperfine splittings between the D and D* mesons. 



For light mesons j25], it has been found that the effects of the bubbles can be large for the pseudoscalar and scalar 
mesons where the additional physics is the anomaly and the glueball respectively, but not for other channels. 
Essentially, the disconnected loops are large if there is additional interesting physics such as glueballs, the anomaly, 
or instantons. It is possible to explore this non-perturbatively from the lattice and we discuss these mechanisms in 
section Ivl 

It is interesting to compare the hyperfine splitting for D mesons with that in charmonium, as there is no contribution 
from the bubble diagrams for the D mesons. In tabled we have collected some results for the D* — D mass splitting 
from quenched QCD. The agreement between experiment and lattice is pretty good for the mass splitting between 
the D* and D. The differences could be explained by the ambiguity in determining the lattice spacing in quenched 
QCD. As noted recently by di Pierro et al. [23|, the hyperfine splitting in the D system is first order in the clover 
coefficient, but the hyperfine splitting in charmonium is second order in the clover coefficient. Hence, the hyperfine 
splitting in charmonium may be more sensitive to the clover coefficient. The differences in the hyperfine splitting may 
potentially be due to remaining errors in the determination of csw- 

It is clearly not sufficient to just assume that the OZI-violating disconnected contributions to charmonium states are 
negligible, particularly as they are responsible for the decay width of rjc of some tens of MeV. In order to clarify these 
issues, we explore from first principles the importance of disconnected contributions to the charmonium hyperfine 
splitting. 



II. SINGLET CORRELATORS 



In lattice studies it is possible to measure separately the non-singlet contribution given by connected correlation 
C{t) (see figure E)) and the flavour singlet contribution which has an additional disconnected correlation D{t) (see 
figure inj . Previous lattice studies have been made of the light pseudoscalar mesons [s^, 113, 0, Issl IsgII and scalar 
mesons [s^, IS HE HO, ^3 ■ For a discussion including some results for vector and axial mesons, see |42j|. 

In the flavour singlet case there is an additional disconnected correlation D{t) to be evaluated. This correlation 
can be written in the form 



D{t)^Nfur,{mL*it)) (1) 

where the disconnected loop 

L{t) = Tr TM-^ (2) 

with the quark propagator and the sum in the trace is over colour, Dirac and spatial indices at time t. Here we 
assume that the hadron under consideration is created by qVq and so the factor of r4 arises from reflection positivity 
(i.e. 74r^ = r4r74). The factor of arises since the Wilson-Dirac fermion matrix Ad is 75 hermitian and hence 
L is real/ imaginary as 75r = r5r75 with = ±1. Since at t = we have that L(0)i*(0) > 0, the disconnected 
correlation -D(O) has sign r^r^. 

At large t where ground state contributions dominate we have 

C{t) = ce~"i* (3) 

and 



C{t) + D(t) = de-"°* 



(4) 
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where mo is the flavour singlet mass and mi the flavour non-singlet mass. Here we are ignoring lighter singlet 
pseudoscalar states (with no charm content, such as rj) which contribute with very small amplitude to C + D. If the 
same meson creation and destruction operators are used for the study of both correlations, with quarks degenerate in 
mass, d and c have the same sign. 

Then by a study of D/C which is given by 

D{t)/C{t) = (d/c)e("i-"'')* - 1 (5) 

the mass splitting between flavour singlet and non-singlet can be explored. Although it might be thought that d = c, 
we have shown previously [i^ that this is not necessarily the case, and indeed sign changes in D/C versus t can be 
required. So, in summary, the slope (increase/decrease) of D/C on a lattice can determine the sign and magnitude 
of mi — toq. For charmonium it is correct to use Nf = 1 in equation ^ since only one flavour of quark can contribute 
to the loops. In our comparisons using lighter quarks for the loops, we will also use Nf = 1. A complete analysis of 
the light two flavour correlator will appear elsewhere [Zlj . 

Since, as we shall see, the disconnected contributions are poorly determined as t increases, it is advantageous to 
remove excited state contributions as far as possible. One technique, pioneered by NefF et al. |43|, is to study the ratio 
of D/C using the ground state contribution to the connected correlator C from a fit. This will be appropriate if the 
disconnected contribution D has only small excited state contributions, as does seem to be the case. 

D{t)/C{t) = (d/c)e('"i-'"«)* - 1 (6) 



III. LATTICE METHODOLOGY 



We use dynamical fermion configurations with Nf = 2 from UKQCD ^j. The sea quarks correspond to k = 0.135 
with a non-perturbative improved clover formalism. The volume was 16'^32. This data set has a scale set by [43 | 
ro/a = 4.754(40) (-1-2 — 90) and pseudoscalar meson to vector meson mass ratio of mp/my = 0.70. Using the value 
To = 0.525 fm then gives — 1790 MeV while the meson mass ratio implies that the sea quarks have masses close 
to that of a strange quark. Wc have already reported the spectrum of the charm-strange mesons and preliminary 
results for the mass of the charm quark on this data set [isl |46{ . 

Local and spatially-fuzzed operators ^47.] are used for meson creation (with two fuzzed links in a spatially symmetric 
orientation with 5 iterations of fuzzing with coefficient given by 2.5*Straight -I- Sum of staples). Thus we evaluate a 
2x2 matrix of local and fuzzed correlators • Mesons created by all independent products of gamma matrices are 
evaluated. 

We measured the connected and disconnected correlations on 201 configurations of size 16'^32 separated by 40 
trajectories for three heavy k values: 0.113, 0.119, 0.125. This data set was used to estimate the k value for the mass 
of the charm quark. Our preliminary estimate for the k value at the charm quark mass was close to 0.119. As the 
aim of this study was to look for the singlet contribution to the charmonium correlators, we did additional runs at 
K = 0.119. At K = 0.119, we computed connected and disconnected correlators separated by 10 sweeps, hence the 
ensemble size was 788. The correlators were then blocked with a block size of 40 sweeps. At k — 0.119, all the results 
reported here are from the higher statistics run. 

For the evaluation of the disconnected correlators, we use 100 stochastic noise sources with the two source trick 
described in |40j . We use sources at every site on the lattice and determine the momentum-zero correlations from 
them. The connected correlators are obtained by explicit inversion from a source (local or fuzzed) 44] . 



IV. ANALYSIS OF THE CONNECTED CORRELATORS 



The lattice spacing used in this data set is large relative to the mass of the charm quark, hence lattice spacing errors 
are a potential concern. In quenched QCD it is computationally possible to use finer lattice spacings, so lattice spacing 
errors can be controlled by "brute force" . The high cost of reducing the lattice spacing in unquenched calculations 
means that a brute force approach will not be feasible for many years with this typ e of fermion action. Hence we 
chose to investigate the heavy quark formalism developed by the Fermilab group [48j. 

The lattice artifacts modify the dispersion [iMIi^ relation: 

= + + 0{p') (7) 

Jyh 

where Mi is known as the "rest mass" and M2 is the kinetic mass (since E = Mi + p'^/(2M2) + 0(p^)). In the 
FNAL lattice heavy quark formalism the rest mass is affected by lattice artifacts, but the M2 mass is the one 



4 



that controls the dynamics of the states. The quahty of the disconnected data precluded us obtaining any useful 
information from the disconnected data with non-zero momentum. A definition of the kinetic mass (M2) in terms of 
the energy E of the meson is: 



M2 dp^ Ip^o 



(8) 



There are a number of different ways to define the quark masses on the lattice. The vector definition of the quark 
mass is 

^v = \{--^) (9) 

where Kcrit is the value of k where the pion mass vanishes. In tree level perturbation theory |48j |. the kinetic definition 
of the the quark mass TO2 is related the vector definition of the quark mass (niy). 

mi = ln(l + mi,) (10) 



1712 171^(2 + 171^) 1 + 

In the ALPHA formulation 50] the vector definition of the quark mass is 0(a) improved using 

Thy = m„(l + hjnrriy) (12) 
where the value of 6^ from perturbation theory is 

6„ = -i- 0.0962g2. (13) 

There are different ways of including the hm term in the calculations. The merits of the different ways are discussed 
in the UKQCD paper {2^ on heavy-light decay constants. The tadpole improved expressions for the FNAL quark 
masses are obtained by replacing ruy with niy/uo, where Ug = l/{8Kcrit)- 

To find the value of k for the charm quark, we interpolate the spin averaged heavy meson mass 

Msav ^ \i3Mv + Mps) (14) 

linearly with the vector definition of the quark mass to the experimental value at 3068 MeV. For the data set with 
correlators separated by 40 sweeps, we fitted a two exponential factorising fit model to the 2 by 2 smearing matrix. 
From that, using Mi, we obtain the kappa value of 0.116 for the charm quark. We comment later on the consequences 
of using other definitions of the mass. 

We used "factorising fits" with three exponentials to fit the two by two matrix of smeared correlators for the higher 
statistics data set. The effective mass plots for the pseudoscalar and scalar channels are in figures and [3 The 
dispersion relation of the heavy-heavy pseudoscalar channel is plotted in figure |21 At the lattice spacing for this 
calculation, the kinetic and rest masses differ by a significant amount. In figure^ we plot the various definitions of 
the meson mass as a function of the quark mass (defined from the vector Ward identity, eq. ISJ. In quenched QCD 
at finer lattice spacings (0.07 fm) UKQCD (2^ has shown that the Mi and M2 masses essentially agree. The high 
computational cost of reducing the lattice spacing forces us to remain in a region where Mi and M2 still differ and 
we determine both masses. We fit equation [7| to the dispersion relation to calculate the M2 mass. The results for the 
various masses are in table Ull 

The correlators with non-zero momentum are noisier than those with zero momentum. There are perturbative 
expressions for the kinetic meson mass in terms of the rest mass of the meson mass. At tree level 

M2^^ = Mi-|-(TO2"mi). (15) 

Our results for the Mi, M2, and Afj^"^ are plotted against the vector definition of the quark mass in figure 01 The 
tree level perturbative expression does not fit the numerical data very well. 

Table Hill contains the mass splittings from the connected correlators. The results for the M2 mass for the excited 
states are very peculiar: we do see a linear behaviour of versus so that M2 can be extracted, but for the excited 



5 



Particle 


reg 


ion 




aMi 


aM2 


My, 


0-+ 


3 - 


13 


2.7/24 


1.549(1) 


2.01(2) 


0.772(8) 


0-+' 


3 - 


13 


2.7/24 


2.02(3) 


1.6(2) 


1.3(2) 


1— 


3 - 


13 


2.3/24 


1.593(2) 


2.06(2) 


0.772(8) 


1—' 


3 - 


13 


2.3/24 


2.09(3) 


1.6(2) 


1.3(2) 


0++ 


3 - 


14 


20/27 


1.790(6) 


2.5(2) 


0.71(5) 


0++' 


3 - 


14 


20/27 


2.38(5) 


- 


- 


1+- 


3 - 


13 


8.8/24 


1.805(9) 






1+-' 


3 - 


13 


8.8/24 


2.31(5) 






1++ 


3 - 


13 


16.4/24 


1.816(7) 








3 - 


13 


16.4/24 


2.39(5) 







TABLE 11: Results for the masses from the fits to the connected data at k = 0.119. The masses of the ground and first excited 
state are shown. The M\ and M2 masses are the meson masses from the Fermilab formalism (from equation|7||. The fit regions 
and /dof are from the momentum zero fits. 
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TABLE IIL Mass splittings from this data set. Sav is the spin averaged mass (see equation 114^ . The experimental numbers 
come from the particle data group. 

states it is smaller than for the ground state. This illustrates the limitation of inter pre ting; M2 as the meson mass 
and implies that a more sophisticated treatment is needed to deal with heavy quarks |5j| . The JLQCD collaboration 
also argue that the kinetic mass is not necessarily superior to the pole mass l5ll |52ll. 

There is a lot of experimental interest in the mass of the r]c{2S) meson Until recently the mass of the 

77c (25) determined from the Crystal Ball collaboration [s^l was larger than the predictions from potential models [55l| . 
However, the new results for the mass of the rjc{2S) from CLEO, BABAR, and BELLE are in much better agreement 
with potential models Is^. In t able II V I we collect the results from some quenched lattice calculations and some recent 
experiments. Our result for the first excited hyperfine splittings is probably effected by lattice spacing errors. We 
discuss the effect of glueballs on the ric{2S) meson in section Ivl 

A. Stochastic noise compared to signal 

We measure the zero momentum disconnected loop L{t) on each time-slice for each gauge configuration. This 
ensemble, for each choice of operator T gives us the values of the standard deviation cTobs given in Table IVl We also, 



Group 


Method 


M^(2S) - M^e(2S) MeV 


Columbia |12] 
CP-PACS fl3] 


anisotropic, lattice 
anisotropic, lattice 


75(44) 
26(17) 


PDG [75]/Belle [56j 
CLEO 
BABAR f58] 
Crystal Ball [54] 


Particle data table 
Experiment 
Experiment 
Experiment 


32 ± 6 ± 8 
43 ± 4 ± 4 
55 ±4 
92 ±5 



TABLE IV: Collection of results for the excited hyperfine splittings between the ip{2S) and rjc{2S) from lattice QCD and 
experiment. 
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TABLE V: Mesons produced by different operators tjiTtp. The standard deviation of the loop operator of eq. |5|is presented. 
Here (Tstoch is the error estimated from the 100 stochastic samples used and this is the used to deconvolute the observed spread 
to give the true standard deviation of the loop (crgaugc). 



^ r t/ f^oha ^stoch f^gauge 

0.135 75 0"+ 33.6 13.91 30.6 

0.135 7fc 1— 14.7 14.45 2.7 

0.135 I 0++ 53.0 15.0 50.8 

0.119 75 0"+ 15.9 8.3 13.6 

0.119 7fc 1— 9.1 9.003 1.2 

0.119 / 0++ 23.6 10.9 20.9 



from our 100 stochastic samples in each case, have the estimate of the standard deviation (Tstoch on the mean of these 
100 samples coming from the stochastic method. We can then deduce the true standard deviation of the gauge time 
slices from (Tgauge = (^obs ~ '''stoch) ■ This is presented in Table Ivl Here the normalisation is such that M = 1 + k 

In an ideal world we would have dstoch ^ '''gauge which would imply that no appreciable error arose from the 
stochastic methods employed. The signal in the vector channel is dominated by the stochastic error. So for that case, 
either many more samples are required or an improved algorithm, such as variance reduction [40l l59j | . For the other 
cases, the stochastic noise is smaller than the intrinsic gauge fluctuation, and so more stochastic samples would not 
improve the results significantly. 

From table the error on the heavy-heavy data is much less than that from the light-light data. This is a 
consequence of increased diagonal dominance of the fermion operator as the mass of the quark is increased. This is 
the basis of improved variance methods ,4^ ,5^ . 

B. Results 

We present in figs. [7| and O some of our results for the ratio of the disconnected correlator to the connected 
correlator for the heavy kappa value and light kappa value respectively. We also measured with fuzzed operators for 
the disconnected diagrams, but they are more noisy than the local operators we present here. 

The error on the disconnected correlator is much larger than that on the connected one. This arises essentially 
because the absolute error on the disconnected correlator stays of the same magnitude as t increases, much as is the 
case for correlations between Wilson loops as used in glueball studies. The connected correlator, in contrast, has an 
approximately constant relative error as t increases. We are forced to consider the ratio of disconnected to connected 
correlator at rather low t-values because of the increasing errors on the disconnected correlator. 

For light (close to the strange quark mass) disconnected contributions we find similar results in fig. |51 to those 
obtained from the lattice previously ,25,] , namely very little signal for vector mesons but a large signal for pseudoscalar 
mesons that will increase the singlet mass over the non-singlet. 

For charm quarks, we present our results in figures [7| and |S1 There is only a statistically significant signal at small 
time values, so no definitive statement can be made. For the pseudoscalar, the slope for non-zero t is positive which 
corresponds to a reduction in the singlet mass compared to the non-singlet. We show on the figure lines corresponding 
to a mass shift of 18 and 36 MeV. This indicates that we cannot rule out a downward shift of the rjc mass by as much 
as 36 MeV. For the vector case, the signal is smaller (in fig. Ej) and shows no sign of a significant slope. 

As we discuss below, we expect the splittings in the vector channel to be small and our results at i 7^ are consistent 
with that. We do find room, however, for splittings of the order of 20 MeV, particularly for the pseudoscalar channel. 
We have shown that the contribution of the singlet correlators to the hyperfine splitting in charmonium may be 
significant. From this calculation, it seems reasonable that the singlet correlators could contribute as much as 20 
MeV to the hyperfine splitting. A more definitive estimate requires more uncorrelated gauge configurations and/or 
improved lattice formulations. 

V. DISCUSSION 

We first recall our previous results from studying the light singlet mesons (2^ on a lattice. The splitting in mass of 
flavour singlet and non-singlet mesons with the same quark content arises from giuonic interactions. The assumption 
that these are small is known as the OZI rule. For the pseudoscalar mesons this splitting is not small (it is related 
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to the 77, T]' mass difference), basically because of the impact of the anomaly. For scalar mesons the splitting is also 
expected to be large because of mixing with the nearby scalar glueball. It is usually assumed that the OZI rule is in 
good shape for the vector and axial mesons and we found small contributions only. 

The picture from the light singlet mesons is that the contribution of the disconnected piece to the correlators is 
small unless there is additional interesting dynamics. In the charmonium system one possible source of the interesting 
dynamics is glueballs. The simplest model is of a flavour singlet state obtained from the mixing of the flavour non- 
singlet state with a glueball, which causes the states to repel in energy, often called an avoided level crossing. We 
would expect this mixing to be strongest when a glueball lay near in energy to the charmonium state and we now 
discuss this. 

Morningstar and Peardon [gJ have computed the excited glueball spectrum in quenched QCD. They obtained 
masses of 2590(40) (130) MeV and 3640(60) (180) MeV for the ground and first excited states of the 0"+ gluebaU 
respectively. Morningstar and Peardon computed the mass of the 1 glueball to be 3850(50)(190) MeV. So it is not 
inconceivable that the r\c mass (2980 MeV) is effected more by glueball states than the J/V' state. In figure^]we plot 
the masses of the glueballs from quenched QCD versus the experimental numbers. With these glueball masses, the 
mixing model predicts that the singlet contribution to the 77c will increase the mass, but that the singlet contribution 
to J will decrease the mass. This would not help to resolve the discrepancy of the charmonium hyperfine splitting 
from non-singlet lattice studies. Moreover this glueball mixing model gives the opposite sign to the pseudoscalar 
mass shift than that indicated by our lattice determination of the disconnected contribution. Bali has also recently 
reviewed the influence of glueball states on the charmonium spectrum _24J . 

The hyperfine splitting between 77c(2S) and tp{2S) states will also be interesting as it may be affecte d by the glueball 
states. The closeness of the glueball state to the ipi^^) state has been noticed by model builders [63 ■ The model 
for hadron decays for vector charmonium states involves the emission of three gluons. This model predicts that the 
branching ratio for '0(25') —^ pir is much larger than experiment. Attempts have been made to use the vector glueball 
and ip{2S) mixing to account for this. If the glueball states have large widths, however, then it is unclear what the 
effect of the states will be on the charmonium spectrum. 

As well as mixing with glueballs, there are other theoretical models which may give guidance on favour singlet 
mass splittings. Isgur and Thackcr discuss the origin of the OZI rule from a quark model and the large Nc limit of 
QCD |42|. Schafer and Shuryak discuss the OZI rule using instanton-based methods j63 |. 

Another approach is to relate the mass splitting to the fact that the decay products (or strongly coupled many-body 
channels) of the singlet and non-singlet state are different. One idea is that a mass shift can arise from the energy 
dependence of the decay width and will be more significant for wider resonances. It will also be possible that mass 
shifts can arise from the back-reaction of the decay channels to the effective propagator. One consequence of this, 
as has been known for a very long time |64j . is that the pole in the complex plane corresponding to a resonance has 
an energy whose real part is lower than the quoted value which corresponds to a phase shift of 90*^. This effect of 
hadron decay on the mass is also an issue for quark model calculations. Isgur and Geiger j6^|6i| discuss the effect 
of decay thresholds on the masses obtained from quark model calculations. In principle, the effect of hadronic decays 
can be introduced into quark models using coupled channel techniques '67*1. However, it is difficult to write down a 
reasonable operator for pair creation. 

This issue of the effect of coupled channels (including open decays) to the mass of a state is one that can be 
illuminated from lattice gauge theory. This is especially so for singlet states, since the lattice enables one to determine 
the relative contributions from the connected and disconnected diagrams separately. 

In the quenched approximation, decays are not treated correctly. This can be a serious problem: the connected 
correlators include only part of the allowed two-body intermediate states and hence anomalous results can be obtained, 
as for the scalar meson 6^. Here we are using a dynamical quark formalism which is explicitly unitary (at least in 
a world with only N f = 2 flavours of quark degenerate in mass). Within this formalism we can add charm quarks 
without expecting any significant breakdown of unitarity from neglected charm quark loops in the vacuum. Then 
the correct treatment of charmonium states is to add the connected and disconnected contributions, as we have 
emphasised. The connected diagram, once one remembers that light quark loops are present in the vacuum, contains 
intermediate states such as DD and D*D etc. It does not contain charmless intermediate states. The hadronic 
decays of those charmonium states below the DD threshold are necessarily to charmless intermediate states. These 
are just the charmless states that are allowed as intermediate states in the disconnected diagram evaluated on the 
lattice. Thus there is a link between the disconnected diagram and the hadronic decay of the charmonium state. For 
an OZI-violating decay, the charm quark and anti-quark must annihilate which is similar process to the contribution 
of the disconnected diagrams (figure O to the singlet correlator. This link is not unambiguous for light quarks: for 
example the substantial 77 -tt splitting (in a world with Nf = 2) arises from the disconnected diagram but no hadronic 
decay of the 77 is allowed energetically. There does not seem to be a simple quantitative link between the mass shift 
caused by the disconnected loop in the correlator and the width of the state. 

This link is explicit in a perturbative treatment of charmonium (even more so for hh): one can evaluate the OZI 
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violating contributions to charmonium from multiple gluon exchange. The pseudoscalar meson allows two gluon 
exchange and so should have much larger effects than for the vector meson where three gluons are needed. Moreover 
the hadronic decays are to multi light-quark states created from these two (or three) gluon intermediate states. 

The computation of strong decay widths from lattice QCD is a hard problem. There are formalisms available, 
but the numerical calculations are quite difhcult. Some of the issues about decay widths and lattice QCD have been 
recently been reviewed [b^ [xOf ■ In a large lattice spatial volume, the effect of coupled two-body decay channels on 
the mass of a state is already taken account of by the formulation, provided one uses a unitary theory with the same 
valence quarks as sea quarks. Thus one should not expect any shift from decay channels. One example of this is that 
the baryon decuplet shows experimentally an equal mass spacing arising from the number of strange quarks present, 
even though the widths of the members vary from 120 MeV (A) to stable (fi). 

On a lattice, at smaller volumes, the two-body momentum states become discrete and this induces small shifts in 
mass. These have been exploited by Liischer to yield information about two-body scatte ring from the lattice. An 
example of a shift in the p mass on a lattice from its coupling to tttt has also been studied [71| ■ 

For hadrons containing only light quarks it is difficult to compute decay widths from first principles on a lattice. 
For charmonium, decay widths can in principle be computed using the NRQCD factorisation formalism "t^, '73^ , or 
from older techniques based on factorising the decay width into a perturbative part and the wave function at the 
origin [741]. 

Thus we should interpret our results as giving an indication of the strength and sign of OZI violating contributions 
to the heavy meson spectrum. These need not correspond to those observed experimentally because we would need to 
extrapolate our lattice results to the continuum limit and to more realistic quark masses (including a third flavour). 
This extrapolation in quark mass could be quite delicate, because of issues such as mixing and decays, as discussed 
above. We note that be mesons will not have these singlet contributions and so the hyperfine splitting for them should 
agree with a lattice calculation using only connected contributions. This may be a useful experimental source of input 
into the composition of such states. 

Both the r^c and J/ip are below the threshold for DD decays. However OZTviolating hadronic decays are allowed. 
The current summary of the rjc properties in the particle data table "tF], quotes the width of the rjc as IGlj'^ MeV. 
However, the latest results for the width of the rjc are larger than the number in the particle data table. CLEO [t^, 
EES 77], and BaBarJp obtain for the ijc width: 27 ± 5.8 ± 1.4 MeV, and 17 ± 3.7 ± 7.4 MeV, and 33 ± 2.5 MeV 
respectively. CLEO '761] note that a larger width (~ 28 MeV) for r/c improves agreement with experiment for the next 
to leading order perturbative QCD expressions |74j for the ratio of the decay width to the two photon width (which 
is more precisely known). This agreement with perturbative estimates also suggests that mixing with glueballs is not 
the dominant mechanism for the rjc decay, and hence for the rjc mass contribution from disconnected diagrams. We 
note that an rjc width of 30 MeV is comparable to the width of some cc mesons above the DD threshold, such as the 
-0(3770) with a width of 24 MeV, although typically the widths of cc mesons above threshold are above 80 MeV. 

The width of the J / tp is accurately known and is 87 keV [t^ . The reason for the smaller width of the vector mesons 
is that the leading order perturbative corrections are 0(02) foi" vector channel, but 0{a2) for the pseudoscalar 
channel. 

Although a width of 30 MeV is small relative to the mass of the rjc, this width is not small relative to the hyperfine 
splitting. It is the mass splittings that are the significant quantities in charmonium. So it is not unreasonable that 
there is a shift coming from OZI-violating intermediate states in the mass of the rjc of the order of 20 MeV and that 
the mass of the J/tp is unaffected. This can be substantiated by more accurate lattice evaluations, our calculation 
leaves room for an effect of such a magnitude but with large errors. 

We have shown that it is necessary to compute the disconnected contributions to the charm correlators to obtain 
accuracies under 3 MeV for mass splittings Q. Studies with anisotropic lattices may be useful to sample more 
intermediate points, but a long reach in physical time is also required. The study of disconnected charm quark loops 
may also be useful for looking at hidden charmonium j79j. 
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FIG. 1: Effective mass plot for the pseudoscalar at k=0.119 




FIG. 3: Dispersion relation for the pseudoscalar at «:=0.119 
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FIG. 4: Different definitions of tlie pseudoscalar meson mass versus tfie vector quark mass, in lattice units. This sliows tfiat 
the central quark mass value (from k = 0.119) is close to the experimental mass using Mi. The perturbative expressions are 
from equation 1151 
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FIG. 5: Connected loop contribution to the propagator 
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FIG. 6: Disconnected loop contribution to the propagator 




FIG. 7: The ratio of disconnected to connected contributions as given by eq. Elat the heavy k value of 0.119. The straight 
lines with slopes of 0.01 and 0.02 (corresponding to a singlet mass 18 and 36 MeV, respectively, lighter than the non-singlet in 
physical units) are drawn to guide the eye. 
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FIG. 8: The ratio of disconnected to connected contributions as given by eq. |n]at the heavy k, value of 0.119. The SESAM 
method is used. The straight lines with slopes of 0.01 and 0.02 (18 and 36 MeV in physical units) are drawn to guide the eye. 




FIG. 9: The ratio of disconnected to connected contributions as given by eq. |^at the light k value of 0.135. 
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S-wave charmonium and glueballs 




FIG. 10: Masses of the experimental pseudoscalar and vector states in charmonium with the glueball masses from quenched 
QCD. 



